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Abstract. We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed 
by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with 
nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer 
(of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic 
weight A > 100) on the top (Haensel & Zdunik, 1990, 2003). The relation between the internal and surface stellar 
temperatures is obtained and fitted by simple expressions. The quiescent luminosity of the hottest (low-mass) and 
coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable 
thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the 
coldest (SAX J1808.4-3658) and the hottest (Aql X-l) neutron stars. The observations of SAX J1808.4-3658 in 
a quiescent state on March 24, 2001 (Campana et al., 2002) can be explained only if this SXT contains a massive 
neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. 
Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal 
emission of SAX J1808. 4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers 
and can serve for testing such models. 
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X-rays: binaries 



1. Introduction 

We study the thermal structure of accreting neutron stars 
in soft X-ray transients (SXTs) — close binaries with a 
low- mass companion (e.g., Chen et al. 1997). Active states 
of SXTs are associated with intense accretion energy re- 
lease and accretion outbursts on the neutron-star surfaces. 
These states are separated by long periods of quiescence, 
when the accretion is switched off or strongly suppressed. 
As noticed by Brown et al. (1998), the spectrum of quies- 
cent emission is well fitted by a neutron-star atmosphere 
model and may thus be of thermal origin, being supported 
by the deep crustal heating due to nuclear transformations 
in the accreted matter. 

Recently the thermal structure and thermal emission 
of neutron stars in the SXTs has been studied by Brown 
et al. (2002), taking into account that hydrogen burning 
in the surface layers may proceed far beyond Fe, up to 
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Te (with nuclear mass numbers A ~ 100), via the rapid 
proton capture process (Schatz et al., 2001). The ashes of 
this burning have large nuclear charges Z, which greatly 
reduces the thermal conductivity of accreted matter and 
increases the internal stellar temperature T in for a given ef- 
fective surface temperature T e g (or a given surface thermal 
luminosity £ 7 ). However, the heavy nuclei may photodis- 
integrate in "superbursts" (Schatz et al., 2003), produc- 
ing nuclei of the iron group, with smaller Z. Moreover, re- 
cently Woosley et al. (2003) have performed new modeling 
of X-ray bursts with updated physics input. Among many 
simulated X-ray bursts, only one anomalous burst pro- 
duced heavy nuclei (A ~ 100), while the other bursts pro- 
duced nuclei with A ~ 60. Taking into account a very wide 
range of physical conditions in bursting neutron stars, we 
consider both possibilities of burning to the elements with 
A ~ 60 and A ~ 100. 

The nuclear ashes, left after bursts and superbursts or 
after steady-state thermonuclear burning in the outermost 
layers, sink in the neutron star crust under the weight of 
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newly accreted matter. With increasing pressure, the sink- 
ing matter undergoes a sequence of nuclear transforma- 
tions (particularly, pycnonuclear reactions), accompanied 
by heat deposition (the so called deep crustal heating in 
accreting neutron stars). Hacnsel & Zdunik (1990) (here- 
after HZ90) studied these processes, starting from the 56 Fe 
ashes (see Bisnovatyi-Kogan 2001, for references to some 
earlier work). Recently the case of ashes of much heavier 
elements has been considered by Haensel & Zdunik (2003) 
(hereafter HZ03), with special attention to 106 Pd ashes. 
The initial mass number A strongly affects the composi- 
tion of accreted matter at densities p < 10 12 g cm" 3 , less 
strongly at higher p, where pycnonuclear reactions oper- 
ate, and moderately affects total crustal heat release (1.45 
MeV and 1.12 MeV per accreted nucleon for A = 56 for 
106, respectively). 

Thermal states of neutron stars in SXTs have been 
studied and compared with observations by a number of 
authors (e.g., Ushomirsky & Rutledge, 2001; Colpi et al., 
2001; Rutledge et al., 2002b; Brown et al., 2002; Yakovlev 
et al., 2003). 

Particularly, Yakovlev et al. (2003) used a simplified 
model of neutron-star thermal structure with the iron 
heat-blanketing envelope, and employed the relation be- 
tween Ti n and L 1 from Potekhin et al. (1997) (hereafter 
PCY). However, even a small mass AM of non-burned ac- 
creted H or He on the surface (AM > (10~ 18 - 10~ 16 ) M 
depending on the surface temperature) noticeably in- 
creases the thermal conductivity of the envelope (Chabrier 
et al. 1997; PCY). On the contrary, the thermonuclear 
burning to high-Z elements (see above) decreases the ther- 
mal conductivity. In the present paper, we take into ac- 
count both effects. First, we derive the relation between 
T; n and L 7 for an envelope composed of a helium layer 
of arbitrary thickness and an underlying heavy-element 
crust, described either by HZ90 or by HZ03 model. Second, 
we use this relation for calculating the thermal states (L 7 
as a function of M) of accreting neutron stars for sev- 
eral neutron-star models with our fully relativistic code 
of neutron-star thermal evolution. We consider five model 
equations of state (EOSs) of matter in the neutron star 
cores for two compositions of this matter — nucleon mat- 
ter and nucleon-hyperon matter. Finally, we compare the 
theoretical results with observations of several SXTs in 
quiescence. We make special emphasis on the hottest and 
coldest neutron stars, in Aql X-l and SAX J1808. 4-3658, 
respectively. 

2. Thermal structure of heat-blanketing envelope 

We calculate the thermal structure of a transiently ac- 
creting neutron star in a quiescent state using the same 
framework as for an isolated neutron star with an accreted 
envelope considered in PCY. A huge energy released in 
the surface layers during the active states of the system 
(heating of the surface by infall of accreted matter and by 
thermonuclear burning of this matter) is quickly carried 
to the surface by the thermal conduction and radiated 



away by the surface photon emission, especially after a 
quiescence onset. Thus, we solve numerically the thermal 
structure equation in the stationary plane-parallel approx- 
imation assuming no energy sources in the heat-blanketing 
envelope. In this way we calculate the temperature distri- 
bution in the blanketing envelope, with the thermal flux 
emergent from warm stellar interiors. The equation to be 
solved is 

dlogT _ 3 PKT C \ 

d logP 16 g T 4 ' [ ' 

where T is the temperature, P is the pressure, g is 
the surface gravity, and K is the opacity. The grav- 
ity is determined by the total gravitational stellar 
mass M and the circumferential stellar radius R: g = 
GM/(R 2 ^Jl-r g /R), where r g = 2GM/c 2 = 2.95 M/M 
km is the Schwarzschild radius. The total opacity K is 
determined by the radiative (photon) opacity K 1 and the 
electron thermal opacity K e : K^ 1 = K~ x + K~ x . The 
effective temperature T c ff is defined by the Stefan law, 

L 7 = 4nR 2 aT* s , (2) 

where a is the Stefan-Boltzmann constant, and L 7 is the 
thermal surface luminosity in a local neutron-star refer- 
ence frame. The apparent luminosity measured by a dis- 
tant observer is = (1 — r g /R) i 7 , and the apparent 
surface temperature inferred by the observer from the ra- 
diation spectrum is T£§ = T e ff a/1 - r g /R (e.g., Thorne 
1977). 

We adopt the standard outer boundary condition to 
Eq. (1) by equating T c s to the temperature T s at the stel- 
lar "radiative surface," which is found from the equation 
P = 2g/(3K). Shibanov et al. (1998) have checked that 
the replacement of this boundary condition by a more re- 
alistic condition, which involves a solution to the radiative 
transfer problem in the stellar atmosphere, has almost no 
effect on the temperature distribution within the blan- 
keting envelope. We assume that the blanketing envelope 
extends to the density pb = 10 11 g cm~ 3 and we integrate 
Eq. (1) from the surface to pb- Thus, we define the internal 
neutron-star temperature as T; n = T(pb). 

We assume further that the outermost neutron-star 
layer is composed of 4 He, which may be left as a non- 
consumed fuel, or may accrete on the stellar surface after 
the last burst. Actually, the outermost layer may be partly 
composed of hydrogen: the replacement of He with H has 
very little effect on the thermal structure of the envelope 
(PCY). We measure the thickness of this light-element 
layer by the parameter 

V = gj 4 AM/M 7 (3) 

where AM is the total mass of He (and H), and g 14 = 
<?/(10 14 cm 2 s _1 ) (<7i4 ranges from <~ 1 to <~ 4, for dif- 
ferent neutron-star masses and EOSs; 514 = 2.43 for a 
"canonical" neutron star of M = 1.4 M Q and R = 10 km). 
The parameter r\ is directly related to the pressure at the 
bottom of the light-element envelope (see, e.g., Ventura 
& Potekhin 2001, Eq. (8)): T H o,ma X = 1.193 x 10 34 r/ dyn 
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cm 2 . Brown et al. (2002) parameterized the thickness of 



the helium layer using the column density 
to r\ as 



it is related 



y = 1.583 x 10 



211 



M 



V 



Mr. " 2 * 2 



(4) 



where Re = -R/10 6 cm. At r\ w 4 x 10 -8 , the helium layer 
reaches the density p « 10 9 g cm" 3 , where it efficiently 
transforms into heavier elements via pycnonuclear reac- 
tions and /3-captures (e.g., Miiller & Langanke 1994). We 
do not consider higher values of r\. 

Under the helium layer, the crust is composed of 
the burst ashes transformed under the action of beta 
captures, emission and absorption of neutrons and (at 
P <L 10 12 g cm~ 3 ) pycnonuclear reactions. We consider 
two model compositions of these layers of the crust: the 
traditional composition originating from the Fe nuclei 
(A — 56, HZ90) and the one originating from the Pd nuclei 
(A = 106, HZ03). 

As in PCY, we employ the OPAL radiative opacity 
tables 1 (Iglesias & Rogers, 1996), interpolated or extrap- 
olated whenever necessary. For the electron opacity (in- 
versely proportional to the electron thermal conductiv- 
ity), we use our code 2 based on the theory presented by 
Potekhin et al. (1999). In order to determine the density 
at a given pressure, we use the Saumon-Chabrier EOS 
for the outer helium layer (Saumon et al., 1995), and the 
EOS of ideal relativistic electron plasma in the deep layers 
of the envelope (using the fitting formulae of Chabrier & 
Potekhin 1998). 

Note that our treatment of the thermal conductiv- 
ity of He and Fe matter takes into account the effects 
of partial ionization, while the thermal conductivity of 
Pd matter (HZ03) is calculated assuming full ionization 
(because of the absence of appropriate ionization mod- 
els). Our estimates show that the partial recombination 
of heavy ion plasma (Z > 30) occurs (very roughly) at 
p < 10 5 g cm- 3 and T < 10 7 K (see Fig. 1 of PCY). If 
the main temperature gradient in the blanketing envelope 
takes place within the indicated domain, our results for 
the Pd/He envelope become inaccurate. We expect, how- 
ever, that the thermal structure of such envelopes will not 
be too different from the structure of equivalent He/Fe 
envelopes (calculated with account for the partial ion- 
ization). Nevertheless, strictly speaking, our calculations 
fail for cold neutron stars (T c g < 1.5 X 10 6 K) with the 
HZ03 envelopes and small amount of He on the surface 
iv ^ 10~ 14 ). We have checked this assumption by con- 
structing artificial models of Hc/Fc envelopes with fully 
ionized Fe. In a cold star with T- ln ~ 6 x 10 6 K and 
AM/M ~ 10~ 20 such a model gives the thermal surface 
luminosity about 2-3 times smaller than the accurate cal- 
culation. We will see (Sect. 3.5) that the effects of partial 
ionization are important for interpretation of observations. 
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Fig. 1. Temperature profiles (temperature versus density) 
through an accreted heat-blanketing envelope of the canoni- 
cal neutron star (M = 1.4 Mq, R — 10 km). The outer part 
of each model envelope (from the radiative surface to a cusp 
of a curve) is composed of 4 He. The thickness of the helium 
layer is measured by parameter r\ — gl^AM/M: r\ = 10 -16 , 
10" 12 , and 10" 8 (marked near the curves). The composition of 
the inner part of the envelope (after the cusp) is taken from 
Haensel & Zdunik (1990). 



Taking each model of the heavy-element accreted en- 
velope (HZ90 and HZ03), we have calculated T- ln for a 
representative sample of ~ 100 pairs of T e g and r\ values, 
in order to find the T e s~T[ n relation. The calculated tem- 
perature profiles are exemplified in Fig. 1 (for i) = 10~ 16 , 
10~ 12 , and 10~ 8 ). As explained in PCY, the light-element 
layer of the envelope (He in our case) provides a lower ther- 
mal insulation than the heavy-element one. This causes a 
sharper temperature drop from the inner isothermal layers 
to the light/heavy-element interface. The interface pro- 
duces a cusp on each profile in Fig. 1. Outside the inter- 
face, the temperature gradient is smaller than just inside. 
Therefore, the thicker the He layer, the higher the effective 
temperature and the emitted radiation flux, for a given 

T in . 

As expected, the T e g-Ti n relation in the HZ90 case is 
almost the same as for a light-element accreted envelope 
on top of the Fe layer considered by PCY (with small cor- 
rections derived by Potekhin et al. 2003). The replacement 
of H by He has little effect on the T e g-T- ln relation; the ef- 
fect is completely negligible if T- ln exceeds the temperature 
of hydrogen thermonuclear burning (~4 x 10 7 K, Ergma 
1986). However, the replacement of the HZ90 crust by the 
HZ03 crust noticeably affects the T c g — T; n relation, as il- 
lustrated in Fig. 2. In this figure, the solid lines and filled 
circles show the dependence of T c g on T; n for the HZ03 
accreted crust, while the crosses and dashed lines show 
this dependence for the HZ90 crust. The symbols demon- 
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Fig. 2. Effective temperature of the canonical neutron star as a 
function of the internal temperature for the HZ90 composition 
of the envelope (dashed lines) and for the HZ03 composition 
(solid lines). 

Table 1. Parameters of the fit (5) 





HZ90 


HZ03 


Ol 


2.87 


2.67 


02 


0.534 


0.560 


<Z3 


2.203 


2.596 


O4 


1.349 


1.136 


o B 


8.328 


9.838 


oe 


0.01609 


0.00967 


a 7 


0.1378 


0.1590 



strate the results of numerical calculations, while the lines 
correspond to the following fitting expressions: 



T e S,6 


= [7T f 4 uU , 6 + (l- 7 )T n V 6 ] 1/4 , 


(5a) 


-^non,6 


- 9l4 a l J in,9' 


(5b) 


TfaU.fi 


= 6.9 gtf C= t tie- 

1 + a 5 J in,9 


(5c) 


7 


= (l + 120T in , 9 y, p=-a 6 /r 1 a ?. 


(5d) 



Here, T eff)6 = T off /10 6 K, T in , 9 = T in /10 9 K, and the pa- 
rameters di (i — 1, 2, . . . , 7) are given in Table 1; T non , 6 
and Tf u n 6 refer to T e gfi in the absence of a helium layer 
and in the presence of most massive helium layer, respec- 
tively. If 5 x 10 6 K < T in < 10 9 K, the difference between 
the calculated and fitted values of T e g does not exceed 4% 
for 77 values used in the calculations: from to 10 -6 for the 
HZ90 composition and from lCT 16 to 1(T 6 for the HZ03 
composition (although the interval of 77, where our physi- 
cal model is applicable, is smaller, as explained above). 



3. Thermal states of transiently accreting neutron 
stars 

3.1. Observations 

We will calculate the quiescent thermal luminosity of 
neutron stars in SXTs versus the mean mass accretion 
rate M (from 1CP 15 to 1CP 9 M Q yr _1 ) and compare the 
results with observations. Let us outline the observations 
first. 

We take the same observational basis as in Yakovlev et 
al. (2003). It consists of five SXTs whose quiescent ther- 
mal luminosity and mean mass accretion rate are mea- 
sured or constrained (Fig. 3): Aql X-l (Rutledge et al., 
2000, 2002a), Cen X-4 (Chen et al., 1997; Rutledge et 
al., 2001), 4U 1608-52 (Chen et al., 1997; Rutledge et al., 
1999), KS 1731-26 (Wijnands et al., 2002; Rutledge et al., 
2002b), and SAX J1808.4-3658 (Bildsten & Chakrabarty, 
2001; Campana et al., 2002). We adopt the same data (see 
Yakovlev et al. 2003 for details) but analyze more care- 
fully the hottest and coldest sources, Aql X-l and SAX 
J1808. 4-3658. The data are rather uncertain and are plot- 
ted thus as thick crosses. The estimated values of M for 
Cen X-l and KS 1731-26 are likely to be the upper limits, 
as indicated by arrows. 

Note that the quiescent luminosity of SXTs may vary 
from one observation to another. It is not clear if this 
variability is associated with the thermal radiation com- 
ponent or the non-thermal (power-law) one, not related 
to the radiation emergent from neutron star interiors. For 
instance, Rutledge et al. (2002a) report temporal variabil- 
ity of radiation from Aql X-l over five months in quies- 
cence after an outburst in November 2000. They fitted 
the thermal component of X-ray spectra with the hydro- 
gen atmosphere model assuming a constant neutron-star 
radius. In the beginning of the quiescent period the ther- 
mal luminosity was estimated as « 9.3 x 10 33 erg s _1 . 
It decreased to about 5.3 x 10 33 erg s _1 in three months, 
increased to about 6.3 x 10 33 erg s _1 in the next month, 
and stayed constant in the last month. In Fig. 3 the ob- 
servational error box is centered at the minimum value, 
L^ 3 « 5.3 x 10 33 erg s _1 , while the upper value of the er- 
rorbar approximately corresponds to the maximum value 
of Recently, Campana & Stella (2003) have proposed 
another interpretation of the same data. They assume that 
the quiescent thermal luminosity stays constant while the 
temporal variability of the Aql X-l radiation is associated 
with the variable nonthermal radiation component (which 
arises outside the neutron star, for instance, due to in- 
teraction of a pulsar wind with an accretion flow from 
the companion star) . Their value of the quiescent thermal 
luminosity (redshifted for a distant observer) is in good 
agreement with the central value of the observational er- 
rorbar in Fig. 3. 

We will also be interested in SAX J1808.4-3658. 
Campana et al. (2002) report its XMM-Newton observa- 
tion in a quiescent state on March 24, 2001, when it was 
extremely weak. The spectrum of its emission was well 
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fitted by power-law. It is most likely that the radiation 
originates in a neutron-star magnetosphere or in a sur- 
rounding accretion disk and is not related to the surface 
thermal emission. Thus, Campana et al. (2002) detected 
no thermal emission but obtained an upper limit on 
in that quiescent state. The limit depends on the model 
of the thermal radiation (black-body or hydrogen atmo- 
sphere model) and on assumed neutron-star parameters 
(mass, radius, distance) and varies from ~ 7.5 x 10 29 erg 
s" 1 to 2.5 x 10 31 erg s" 1 . 

We consider two possible values (Figs. 3 and 6) of the 
upper limit on discussed by Campana et al. (2002). 
They seem to be more realistic than several other pos- 
sible upper limits mentioned by these authors. The first 
value, = 1.5 x 10 30 erg s _1 , is the upper limit on 
the unabsorbed bolometric thermal luminosity with the 
black-body (BB) spectrum. It is obtained from the 90% 
upper limit on the unabsorbed bolometric black-body flux 
of ~ 2 x 10~ 15 erg s _1 cm~ 2 . To translate the flux to the 
luminosity we adopt the distance to the source, 2.5 kpc, 
cited by these authors. The second value, = 7.5 x 10 29 
erg s _1 , is the upper limit on the unabsorbed bolometric 
thermal luminosity obtained with the neutron-star hydro- 
gen atmosphere (NSA) model. 

As mentioned, e.g., by Campana et al. (2002) and 
Yakovlev et al. (2003), the quiescent thermal luminosity 
of SAX J1808. 4-3658 is so low that the source should con- 
tain a very cold neutron star (with strong neutrino emis- 
sion). Particularly, as shown by Yakovlev et al. (2003), 
the data are compatible only with the model of a massive 
neutron star whose core contains nucleons (and possibly 
hyperons) and does not contain any exotic matter (pion- 
or kaon condensates or quark matter, which would im- 
ply lower neutrino emission and a noticeably hotter star). 
This restricts the class of possible EOSs of dense matter 
to the models of nucleon/hyperon matter. We extend the 
analysis of Yakovlev et al. (2003) by simulating thermal 
states of neutron stars using an exact cooling code and 
several models of nucleon/hyperon EOSs. 



3.2. Numerical simulations 

Following Brown et al. (1998), we assume that neutron 
stars in SXTs are warmed up by deep crustal heating in 
accreted matter. It is supposed that the heat released at 
their surface in active states is radiated away, especially 
during quiescence, and does not warm the stellar interiors. 
Neutron stars in SXTs are thermally inertial objects, with 
typical relaxation times <~ 1 — 10 kyr (Colpi et al., 2001), 
while the mass accretion rate M in SXTs varies on much 
shorter time scales. Therefore, we will study a global ther- 
mal state of a transiently accreting neutron star by replac- 
ing a variable deep crustal heating with the time-averaged 
heating determined by the time-averaged accretion rate 



M = (M). The deep-heating power is 
QM 



^•° 3Xl ° 33 ( lO-i°M yr-i ) GJv) "S^W 

where ton is the nucleon mass, and Q is the total amount 
of heat released per one accreted nucleon (Q = 1.45 MeV 
and 1.12 MeV for the HZ90 and HZ03 models, respec- 
tively). 

We simulate thermal states of accreting neutron stars 
with our fully relativistic code of neutron-star thermal 
evolution (Gnedin et al., 2001) by solving the stationary 
thermal-balance equation: 

LS° h (M) = L~(T in ) + L~(T eff ) ) (7) 

where is the neutrino luminosity of the star as de- 
tected by a distant observer, and L°^ h is the deep heating 
power for the same observer (see, e.g., Yakovlev et al. 2003 
for details). Since the mass accretion rates in SXTs are de- 
termined with large uncertainties, we do not make a dif- 
ference between L dh and in Eq. (7). The neutron star 
interior (p > pb) is assumed isothermal (with the proper 
account of the general relativistic effects) . In this case L£° 
is a function of T; n , while T- m is related to T e g (Sect. 2). 

We have updated our code in three respects. First, we 
have incorporated the internal energy sources associated 
with deep crustal heating. Second, we have modified the 
relation between the surface and internal stellar tempera- 
tures in accordance with the results of Sect. 2. Third, we 
have included new microphysics (neutrino emissivities and 
heat capacities) which allows us to consider neutron-star 
cores containing A, S°, E~, £+, S~, and S° hyperons (in 
addition to neutrons, protons, electrons, and muons). 

The code calculates heating curves, the redshifted sur- 
face thermal luminosity of the star (or the effective 
surface temperature T£§) versus the mean mass accretion 
rate M, to be compared with observations. Examples are 
shown in Figs. 3-6. 

3.3. Model equations of state of dense matter 

We will use five model EOSs in neutron star cores. Three 
of them (EOSs Nl, N2, and N3 listed in Table 2) refer to a 
nucleon dense matter (neutrons, protons, and electrons), 
while other two (EOSs NH1 and NH2 listed in Table 3) 
refer to a matter containing nucleons, electrons, muons 
and hyperons (of all types). 

EOSs Nl, N2, and N3 are modifications of the phe- 
nomenological EOS proposed by Prakash et al. (1988). 
All of them allow the direct Urea process (Lattimer et al. , 
1991) to operate in a sufficiently dense neutron-star mat- 
ter, p > po- The threshold density po is given in Table 
2. We present also the parameters of the maximum-mass 
configurations (mass M max , central density p Cl and radius 
R) and the parameters of low-mass neutron-star config- 
urations (with M = 1.1 M© as an example). EOS Nl 
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Table 2. EOSs of nucleon matter; low-mass and maximum-mass neutron star configurations 







Low-mass 


Maximum-mass 


EOS 


PD 


M 


Pc 
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M 


Pc 


R 




(g cm" 3 ) 


(Me) 


(g cm -3 ) 


(km) 


(M e ) 


(g cm" 3 ) 


(km) 


Nl 


7.85 x 10 i4 


1.1 


6.23 x 10 i4 


13.20 


1.977 


2.578 x 10 ib 


10.75 


N2 


1.298 x 10 15 


1.1 


8.50 x 10 14 


12.18 


1.73 


3.25 x 10 15 


9.71 


N3 


1.269 x 10 15 


1.1 


1.217 x 10 15 


11.31 


1.46 


4.04 x 10 15 


8.91 



Table 3. EOSs of hyperon matter; maximum-mass hyperonic stellar configurations 











Maximum-mass 


EOS 


Ph 


Pd 


pDl 


M 


Pc 


R 




(g cm" 3 ) 


(g cm" 3 ) 


(g cm" 3 ) 


(M e ) 


(g cm" 3 ) 


(km) 


NH1 


5.495 x 10 14 


4.06 x 10 14 


>pr x 


1.975 


2.35 x 10 1S 


11.41 


NH2 


5.24 x 10 14 


5.25 x 10 14 


2.2 x 10 15 


1.740 


3.891 x 10 15 


9.089 



implies model I of the symmetry energy and the com- 
pression modulus of saturated nuclear matter K = 240 
MeV. EOSs N2 and N3 imply, respectively, K = 180 MeV 
and K = 120 MeV, and the symmetry energy proposed 
by Page & Applegate (1992). EOS Nl is rather stiff: it 
yields M max « 2M Q . EOS N2 is softer, with smaller M max . 
The symmetry energy is overall smaller, which results in 
a lower proton fraction and higher direct Urea threshold. 
EOS N3 is an example of a quite different, soft EOS, with 
M max = 1.46M . 

EOSs NH1 and NH2 in Table 3 refer to hyperon mat- 
ter. The hyperons appear at sufficiently high densities; the 
density of the appearance of the first hyperon is denoted 
by ph- In the inner cores of massive hyperonic stars, a 
variety of direct Urea processes are open. They involve 
nucleons and hyperons, electrons and muons (Prakash et 
al., 1992) and give the major contribution to the neu- 
trino luminosity of massive stars. In Table 3, pd means 
the threshold density for the first direct Urea process (hy- 
peronic or nucleonic). However, if the density of hyper- 
onic matter grows to essentially supranuclear values, the 
fraction of leptons (electrons and muons) becomes lower 
(they arc replaced by X" and S" hyperons). Since the lep- 
tons are important participants of Urea processes, their re- 
duced fraction may result in switching off these processes 
at high density. In Table 3, pdi is the highest density of 
operation of the last direct Urea process. 

EOS NH1 is given by the relativistic mean field model 
3 of Glendenning (1985) and may be unrealistically stiff 
for hyperonic matter (M max « 2M ). The first hyperon, 
A, appears at p = p^, where nucleon direct Urea processes 
are already on (p D < p h ). The reduction of lepton fraction 
at high densities is not too strong. As a result, the direct 
Urea processes still operate at the center of the maximum- 
mass star. 

EOS NH2 is EOS2 NASS of Balberg et al. (1999). It is 
softer: M max « 1.74 M©. The first hyperon, S~, appears 
at the density ph which is only slightly lower than po- 
The first direct Urea process, which opens with increasing 



p, involves hyperons. A very strong reduction of lepton 
fraction at high densities switches off all direct Urea pro- 
cesses at nearly the same density pdi hi central kernels of 
massive stars. 

3.4. Limiting models 

Realistic models of neutron stars should take into account 
possible superfluidity of baryons (nucleons and hyper- 
ons) in stellar interiors (e.g., Lombardo & Schulze, 2001; 
Balberg et al., 1999). Microscopic calculations give a large 
scatter of critical temperatures of baryon pairing, depend- 
ing on a baryon-baryon interaction model and many-body 
theory employed. All these calculations indicate that su- 
perfluidity should disappear in essentially supranuclear 
matter (in the centers of massive neutron stars). 

The theory predicts the existence of three types of 
transiently accreting neutron stars with different thermal 
structure (e.g., Yakovlev et al. 2003): low-mass, medium- 
mass, and high-mass stars. Low-mass stars have so low 
central densities that the direct Urea processes are not 
operative, being cither strictly forbidden (p c < pu) or to- 
tally suppressed by a superfluidity (which may be strong 
at not too high densities) . Accordingly, these stars have a 
low neutrino luminosity and are hottest, for a given mass 
accretion rate. The central densities of high-mass neutron 
stars are noticeably higher than pd and the density, where 
the superfluidity dies out and ceases to suppress the di- 
rect Urea process. These stars possess inner cores with 
fully open direct Urea processes and greatly enhanced neu- 
trino luminosity. Hence, the massive stars are coldest, for a 
given accretion rate. Medium-mass stars are intermediate 
between the low-mass and high-mass ones. 

Thus, the highest heating curves (Fig. 3) correspond to 
low-mass stars, and the lowest curves to high-mass stars. 
Increasing the stellar mass from the lowest to the highest 
values, one obtains a family of heating curves of medium- 
mass stars which fill in space between the highest and 
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Table 4. Variation limits of three parameters of heating mod- 
els 



Limiting parameter 


Limit 1 


Limit 2 


Mass M 


1.1 Mq 

(low-mass star) 


M max (high- 
mass star) 


Accreted crust model 


HZ90 


HZ03 


Mass of He layer, 
AM/M = V /g 2 14 





itr 8 



lowest curves. Any observational point between the lowest 
and highest heating curves can be explained by employed 
neutron-star models. 

Since the problem involves many parameters, we re- 
strict ourselves to the limiting cases. They give the high- 
est or lowest heating curves and constrain thus theoret- 
ical values of L2°. The limiting low-mass and high-mass 
neutron-star configurations seem to be more robust with 
respect to physics input (EOS and superfluid properties of 
stellar cores) than the medium-mass configurations. This 
has been proven for cooling isolated neutron stars (e.g., 
Kaminker et al. 2002) and has to be true for transiently 
accreting neutron stars with deep crustal heating because 
of the direct correspondence of cooling and heating prob- 
lems (Yakovlev et al., 2003). Specifically, we employ lim- 
iting values of three parameters listed in Table 4. 

The first parameter is the stellar mass; the limiting 
values refer to low-mass and high-mass neutron stars. For 
representative low-mass models, wc take the ones with 
M = 1.1 M Q and nucleon cores (Table 2). We assume 
further the strong proton superfluidity in the cores of low- 
mass stars, with the critical temperature of proton pairing 
equal to 5x 10 9 K. This superfluidity totally suppresses the 
modified Urea process, so that the neutrino luminosity of 
the star is determined by neutron-neutron bremsstrahlung 
process (just as for cooling neutron stars; e.g., Kaminker 
et al. 2002). Actually, our models of low- mass stars are 
fairly insensitive to a specific model of the proton super- 
fluidity as long as the superfluidity is sufficiently strong to 
suppress the modified Urea process. 

Recently the neutrino emission in neutron-neutron 
bremsstrahlung has been reconsidered by van Dalen et 
al. (2003). The authors use several realistic models of 
neutron-neutron interaction and conclude that the emis- 
sivity of the process is about four times lower than in 
the simplified one-pion exchange (OPE) model. Our code 
incorporates the results of Friman & Maxwell (1979) ob- 
tained using the OPE model with phenomenological cor- 
rection factors. We have checked that the emissivity pro- 
vided by the code agrees with the improved results of van 
Dalen et al. (2003), i.e., it is already reduced by a factor 
of ~ 4 with regard to the OPE model of van Dalen et al. 
(2003). 

Furthermore, we assume a weak triplet-state neutron 
pairing in neutron-star cores, with maximum critical tem- 



perature < 2 x 10 8 K. This pairing has actually no ef- 
fect on thermal states of neutron stars and can be ig- 
nored. A stronger neutron pairing would produce a pow- 
erful neutrino emission and fast cooling of isolated neu- 
tron stars (e.g., Kaminker et al., 2002; Yakovlev et al., 
2003b), in sharp disagreement with observations of these 
objects. For simplicity, we neglect the effects associated 
with singlet-state neutron pairing in neutron-star crusts 
(e.g., Yakovlev et al. 2001; Potekhin et al. 2003). To rep- 
resent the high-mass stars, we use the maximum-mass con- 
figurations (Tables 2 and 3). The effects of superfluidity 
are expected to be minor in such a star, and we neglect 
them in our calculations. 

The second and third parameters in Table 4 test the 
composition of accreted envelopes. The second parameter 
specifies the composition of deep accreted crust; we con- 
sider two limiting cases, the HZ90 or HZ03 crusts (Sects. 
1 and 2). The third parameter is the relative mass of the 
surface helium layer (Sect. 2). Our two limiting cases are: 
AM/M = (no He layer) and AM/M = 10~ 8 (most 
massive He layer). We will see that the effect of the he- 
lium layer on saturates with increasing AM/M, ow- 
ing to which the limit of most massive envelope is actually 
achieved at AM/M < 10~ 8 . 

3.5. Overall results (EOS Nl) 

We start with the overall analysis, taking EOS Nl in the 
stellar cores as an example. The effects of different EOSs 
will be analyzed in Sects. 3.6 and 3.7. Theoretical models 
are compared with the observations (Sect. 3.1) in Fig. 3. 
Four upper heating curves refer to a low-mass neutron 
star, while three lower curves refer to a high-mass one. The 
solid and dashed curves are for the HZ03 crust, while the 
dash-and-dot and dotted curves are for the HZ90 crust. 
The solid and dashed-and-dot curves are calculated as- 
suming no He on the surface, while the dashed and dotted 
curves are for the massive He layer. We do not present the 
solid curve (the HZ03 crust with the Pd composition ex- 
tended to the surface) for the maximum-mass star. In this 
star, the effects of partial ionization of Pd plasma would 
be very important but they are neglected in our analysis. 
We expect, however, that this curve is not too different 
from the dash-and-dot one (see below). 

Our calculations confirm the main results of the previ- 
ous studies (e.g., Yakovlev et al. 2003): one can explain the 
data, assuming that neutron star cores are composed of 
nucleons; Aql X-l and 4U 1806-52 may be treated as SXTs 
containing low-mass neutron stars; Cen X-4 and KS 1731- 
26 may be treated as SXTs containing medium-mass neu- 
tron stars; SAX J1808. 4-3658 seems to be a source with 
a high-mass star; all these neutron stars cool mainly via 
neutrino emission from their cores (their photon thermal 
luminosity is much weaker than the neutrino one) except 
possibly Aql X-l, whose regime is intermediate between 
the neutrino and photon cooling. In the neutrino regime, 
Eq. (7) reduces to Lf h (M) w L~(T in ), which yields the 
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Fig. 3. Theoretical quiescent thermal luminosity of neutron 
stars (with EOS Nl in the cores) versus mean mass accre- 
tion rate confronted with observations of five SXTs. The larger 
upper limit of the luminosity of SAX J1808. 4-3658 assumes 
the blackbody (BB) surface emission; the smaller limit - lower 
horizontal bar - is obtained with the neutron-star atmosphere 
(NSA) model. The upper curves refer to the low- mass (1.1 Mq) 
star (with strong proton superfluidity in the core), while the 
lower curves to the maximum-mass star. Four curves for the 
low-mass star are calculated assuming HZ90 or HZ03 accreted 
crusts, without or with a massive helium surface layer; singly 
hatched region shows ranges of variations of quiescent luminos- 
ity due to different mass of helium layer for the HZ03 crust; 
double hatched region is the same for the HZ90 crust. The 
heating curves for the high-mass star are the same but the 
solid curve (for the HZ03 crust without helium layer) is ab- 
sent because of the partial ionization problem (see text); it is 
expected to be close to the dash-and-dot curve. 



internal temperature T- ln , while the surface thermal emis- 
sion is adjusted to this T- m (see Yakovlev et al. 2003 for 
details). 

Our present results enable us to extend the consider- 
ation of Yakovlev et al. (2003). Since we are mainly in- 
terested in the neutrino cooling regime, the internal tem- 
peratures of neutron stars of a given mass and accretion 
rate are determined by neutrino emission and deep crustal 
heating model (HZ90 or HZ03) and do not depend on 
the presence of He on the surface. On the other hand, 
the T e ff — Ti n relation is actually the same for the HZ90 
and HZ03 scenarios, if massive layers of light elements are 
present. A small difference between the dashed and dot- 
ted curves is solely determined by different amount of heat 
released in the HZ90 and HZ03 crusts, see Eq. (6). The 
heat released in the HZ90 crust is slightly larger, so that 
the dotted curves are slightly higher than their dashed 
counterparts. A larger difference between the HZ90 and 
HZ03 crusts without a light-element layer is mainly deter- 



mined by different thermal insulations of heat-blanketing 
envelopes (Sect. 2). The high-Z (HZ03) heat blanketing 
envelope is less heat transparent. Hence, the surface tem- 
perature is smaller and the solid curve is lower than its 
dot-and-dashed counterpart. Still larger difference occurs 
between the scenarios with and without helium layers. A 
massive helium layer is much more heat transparent than 
the layer of HZ90- or HZ03-matter. Accordingly, the dot- 
ted and dashed heating curves go noticeably higher than 
the associated dash-and-dot and solid curves. 

As seen from Fig. 3, the presence of light elements on 
the surface of Aql X-l simplifies theoretical treatment of 
Aql X-l as an SXT containing a low-mass neutron star. 
The spectrum of the object is well described by hydrogen 
atmosphere models, which is in line with the assumption 
that the neutron star has the surface layer of light ele- 
ments. The effects of different EOSs for the interpretation 
of this source are discussed in Sect. 3.6. The neutron star 
in 4U 1608-52 may be treated either as a low-mass neu- 
tron star without a massive light-element layer, or as a 
medium-mass neutron star with such a layer. 

The interpretation of SAX J1808.4-3658 is of special 
interest. EOS Nl adopted in Fig. 3 is consistent only 
with the black-body thermal emission (with no massive He 
layer on the stellar surface) and disagrees with a neutron- 
star-atmosphere spectrum. Although the upper limits of 
inferred with the black-body and neutron star atmo- 
sphere models seem to be not very certain, our results 
indicate that the neutron star in SAX J1808. 4-3658 is so 
cold that it is barely explained by the theory. In Sect. 3.7 
we will show that the theoretical explanation is relaxed if 
the neutron star contains a hyperonic core. In any case, 
the star should have no massive layer of light elements on 
the surface. 

The limiting cases of no He layer and massive He layer 
deserve special comments. The mass of the light-element 
layer may vary from one quiescent stage to another or 
even during one quiescent stage due to residual accretion. 
Therefore, the heating curves for a massive He layer and 
without it represent the upper and lower limits of the qui- 
escent thermal luminosity of the same star (Brown et al., 
2002). 

Figure 4 shows the surface thermal luminosity of neu- 
tron stars with the HZ90 and HZ03 crusts versus the mass 
AM of the helium layer. The upper curves refer to the low- 
mass star with the mean accretion rate of Aql X-l. The 
lower curves correspond to the high-mass star with the 
mean accretion rate of SAX J1808. 4-3658. The growth of 
L 1 ^ with increasing AM saturates, so that the limit of very 
massive He layer is actually achieved at AM <C 10~ 8 Mq. 

Since the neutron star in Aql X-l is hot, the main 
temperature gradient occurs in deep layers of its heat- 
blanketing envelope. One needs at least AM ~ 10~ 13 Mq 
of helium to affect the thermal structure, and AM ~ 
10~ 10 Mq of helium (extending to p ~ 10 7 g cm~ 3 ) to 
achieve the limit of the most massive He layer. The in- 
crease of AM enhances the thermal luminosity by a fac- 
tor of 2.5-3, in agreement with the results by Brown et 
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Fig. 4. Theoretical quiescent thermal luminosity of neutron 
stars (with EOS Nl in the cores and the HZ90 or HZ03 crusts) 
in Aql X-l and SAX J1808.4-3658 (as in Fig. 3) versus the mass 
of the surface He layer. The dotted lines refer to oversimplified 
(inaccurate) models which neglect partial ionization of heavy 
elements (Fe or Pd) in outermost layers. 



al. (2002). It is instructive to note that the mass of the 
neutron-star atmosphere in Aql X-l at the optical depth 
r = 2/3 would be - 10~ 24 M , while at r = 10 it would 
be ~ 1O~ 2O M0. Therefore, the spectrum of thermal ra- 
diation formed in the stellar atmosphere is affected by a 
much smaller amount of light elements than the thermal 
structure of the neutron star envelope. Note also that the 
plasma remains almost fully ionized in all surface layers 
of this hot neutron star. 

The situation with SAX J1808.4-3658 is different. 
The neutron star is much colder, and the main temper- 
ature gradient shifts to the surface. Variations of 
with increasing AM reach one order of magnitude. Even 
AM ~ 10 _20 -Mq of He (comparable with the mass of 
the atmosphere at r ~ 10, less than 1 cm under the 
surface) affects the surface thermal luminosity, while the 
limit of the most massive envelope is nearly achieved at 
AM ~ (10~ 14 - 10~ 13 ) M Q . At these AM heavy elements 
are partially ionized. As explained in Sect. 2, we have 
taken into account the partial ionization of Fe but not 
Pd. Assuming artificially full ionization of Fe or Pd, we 
obtain the dotted curves in Fig. 4 (the dotted curve for 
the Pd/He crust at AM < 10~ 15 Af is plotted from direct 
calculation rather than from the fit expressions, Eqs. (5)). 
This approximation is seen to be rather inaccurate for the 
Fc/He envelope at AM < 10 _14 M Q , and is expected to 
be even more inaccurate at these AM for the Pd/He en- 
velope. In this respect, cold neutron stars (particularly, 
SAX J1808. 4-3658) can serve as laboratories for study- 
ing ionization equilibrium of dense matter, a complicated 
theoretical problem whose solution is model dependent. 



3.6. Low-mass neutron stars: Aql X-l 



The curves in Figs. 3 and 4 are calculated for one model 
EOS in a neutron star core. The calculations with all 
five EOSs listed in Tables 2 and 3 reveal that the heat- 
ing curves of low-mass and high-mass neutron stars are 
not too sensitive to these EOSs. The heating curves of 
medium-mass stars do depend on the EOS (just as for 
isolated cooling neutron stars, e.g., Kaminkcr et al. 2002) 
which will be studied in the next publication. Here, we 
restrict ourselves to the limiting low-mass and high-mass 
models. 

Let us start with the low-mass stars. Figure 5 shows 
the heating curves of 1.1 M neutron stars with the cores 
described by EOSs Nl, N2, and N3 and the envelopes 
composed of either the HZ90 or HZ03 matter, without 
and with He layers. The left panel displays our tradi- 
tional heating curves, L^(M), while the right panel gives 
T°g(M) for the same scenarios. 

The left and right panels of Fig. 5 are seen to be dif- 
ferent. If we fix a model of accreted crust and the value 
of AM/M, the L^°(M) heating curves are almost uni- 
versal - insensitive to the EOS in the stellar core. We 
plot them by wide hatched lines in the left panel: the 
curves for different EOSs would be almost indistinguish- 
able. In addition, we have calculated the heating curves for 
a M = 1.3 M Q neutron star with EOS Nl (not presented 
in Fig. 5) and check that they are also indistinguishable 
from corresponding 1.1 Mq curves. On the other hand, the 
T°£ (M) heating curves noticeably depend on the EOS: the 
luminosity L^(M) is the same, but stellar radii are dif- 
ferent (Table 2), hence different T e <g(M). Notice that the 
theory of cooling isolated neutron stars predicts just the 
opposite property of the cooling curves (L?° or ver- 
sus stellar age t): the cooling curves T^(t) of low-mass 
neutron stars are almost universal (e.g., Kaminker et al. 
2002), whereas the curves L™(t) should noticeably depend 
on the EOS. 

The upper heating curves in Fig. 5 correspond to the 
stars with massive light-element layers, while the lower 
curves are calculated assuming no such layers. The dif- 
ference between the lower and upper curves is quite pro- 
nounced. As seen from Fig. 5, observations of Aql X-l 
with simultaneous stringent determinations of and 
T£g have potential to constrain the EOS in the core of 
low-mass neutron stars and to determine the composition 
of the surface layers. On the other hand, the neutron star 
is so hot that it is nearly at the edge of theoretical ability 
to explain hot stars. Were the quiescent thermal luminos- 
ity > 10 34 erg s^ 1 detected in future observations of 
Aql X-l, it could not be explained by deep crustal heat- 
ing. This may be regarded as an additional test of the deep 
crustal heating mechanism. 
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Fig. 5. Quiescent thermal luminosity (left panel) or effective surface temperature (right panel) vs. M for a 1.1 M© star with 
a nucleon core (containing strongly superfluid protons) confronted with observations of Aql X-l. Set of four heating curves: 
with HZ90 or HZ03 crusts, with and without massive helium layer, is shown for the three EOSs in the stellar cores (EOSs Nl, 
N2 and N3). Fixing the accreted crust model, we obtain almost the same luminosities (the curves for EOSs Nl, N2, and N3 
almost coincide; left panel) but different effective temperatures (right panel). Hatched regions in the left panel display ranges 
of possible variations of due to variable thickness of the He layer for the HZ90 and HZ03 crusts. 



3. 7. The coldest massive neutron star: SAX 
J 1808. 4-3658 

Let us turn to the high-mass neutron-star models. Figure 6 
displays the heating curves of maximum-mass neutron 
stars with five EOSs in their cores, the HZ90 or HZ03 
composition of the crust, with or without helium surface 
layers. If we fix the composition of the crust and the he- 
lium mass, then the difference between the thermal lumi- 
nosities of neutron stars with four EOSs in the core (nu- 
cleon EOSs Nl, N2, and N3, and hyperon EOS NH2) is 
small. Thus, just as for low-mass neutron stars (Sect. 3.6), 
the heating curves (M) are nearly universal (while the 
curves T^(M) would be much more different). 

Employing the four EOSs (Nl, N2, N3, and NH2), we 
come to the same conclusions as in Sect. 3.5: the theory 
is consistent with the upper limit of inferred using 
the blackbody model (without any light-element layer on 
the surface), and it is inconsistent with the limit of 
inferred with the neutron-star atmosphere model. The ap- 
pearance of a helium layer can raise the quiescent thermal 
luminosity up to an order of magnitude. 

Nevertheless, the stars with the hyperon EOS NH1 
disobey the universality: they are noticeably colder. Their 
coldness is explained by very strong neutrino emission 
from from the hyperonic stellar core. As discussed in Sect. 
3.3, the fraction of leptons in the central region of the 
maximum-mass neutron star with EOS NH1 is sufficient 
to keep all direct Urea processes open everywhere within 
the inner core, including the stellar center. The neutrino 
emission becomes exceptionally intense, and the star very 
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Fig. 6. Quiescent thermal luminosity vs. M for maximum- 
mass neutron-star models confronted with observations of SAX 
J1808.4-3658. Based on 15 curves for five EOSs in the stellar 
core in the presence of massive helium surface layers (HZ90 
and HZ03 crusts, hatched lines) and in the absence of these 
layers (HZ90 crust, solid lines). The EOS types are marked to 
the left of the lines. The curves for EOSs N2 and NH2 are al- 
most indistinguishable. Taking the massive helium layer and a 
fixed EOS, we have two almost coincident curves for HZ90 and 
HZ03 crusts, plotted as a wide hatched line. 
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cold. Such a cold high-mass hyperonic star without a mas- 
sive helium surface layer is in much better agreement with 
the observations of SAX J1808. 4-3658. As seen from Fig. 
6, in this case the theory is not only consistent with the 
upper blackbody limit on L^°, but becomes in reasonable 
agreement with the upper limit inferred using the atmo- 
sphere model (which is not the case for other EOSs). This 
opens a potential possibility to discriminate between the 
different nuclcon/hypcron EOSs and gives a tentative in- 
dication that the neutron star in SAX J1808. 4-3658 con- 
tains a hyperonic core with a not too low lepton fraction. 

As discussed in Sect. 3.3, the lepton fraction for EOS 
NH2 decreases very rapidly with increasing density. This 
switches off all direct Urea processes at p > pr>i, in the 
central kernel of the maximum-mass neutron star, reduc- 
ing thus the neutrino emission to the level of neutron stars 
with nucleon cores. 

Thus, further observations of SAX J1808.4-3658 in 
quiescent are highly desirable. Future detections or con- 
straints of the quiescent thermal luminosity would be 
most important, especially those which give lower values 
of than the present ones. For instance the value of 
« 6 x 10 29 erg s _1 would definitely imply a massive 
hyperon neutron star with not too small fraction of lep- 
tons in its center (to keep the direct Urea processes open) 
and without any massive surface layer of light elements. 
Lower values of could not be explained by the current 
model. 

4. Conclusions 

We have considered (Sect. 2) the growth of temperature 
within the heat-blanketing envelope of a transiently ac- 
creting neutron star in a quiescent state. We have ana- 
lyzed two basic models HZ90 and HZ03 of the accreted 
crust, calculated by Haensel & Zdunik (1990, 2003). In 
all cases we consider the possible presence of a thin 
(AM < 10 _8 M Q ) layer of light elements (H or He) on 
the surface. We have calculated the relations between the 
internal and surface temperatures of neutron stars with 
the HZ90 or HZ03 crusts and fitted the results by simple 
expressions. 

Using these results, we have modeled (Sect. 3) ther- 
mal states of transiently accreting neutron stars in SXTs, 
assuming that these states are regulated by deep crustal 
heating in accreted matter. We have considered five model 
EOSs of nucleon or nucleon-hyperon matter in neutron 
star cores, representative models of low-mass and high- 
mass neutron stars, HZ90 and HZ03 models of stellar 
crusts, without light elements on the stellar surfaces or 
with maximum amount of light elements. The results give 
the upper and lower limits of the quiescent thermal lu- 
minosity of SXTs, depending on the amount of light ele- 
ments at the neutron-star surface in a particular quiescent 
period. 

We have compared the theory with observations of five 
SXTs. The most important are two sources, Aql X-l and 
SAX J1808. 4-3658. Aql X-l can be treated as a low-mass, 



warm neutron star. Its future observations may constrain 
the EOS in the nucleon core of a low-mass neutron star, 
elucidate the composition of accreted matter, and test 
the deep crustal heating hypothesis. The second source, 
SAX J1808. 4-3658, can be treated as a very cold mas- 
sive neutron star with nucleon or nucleon-hyperon core; 
a hyperonic core with not too low fraction of electrons is 
more preferable. Future observations of SAX J1808. 4-3658 
in quiescence may enable one to distinguish between the 
EOSs in massive nucleon/hyperon stellar cores and check 
the ionization models of heavy-element plasma in surface 
layers of neutron stars. 

The assumption that the quiescent thermal emission of 
SXTs is produced by the deep crustal heating (Brown ct 
al., 1998) remains still a hypothesis. However, the theory 
of deep crustal heating (Haensel & Zdunik, 1990, 2003) is 
solid: accreting neutron stars should be heated from inside 
and this effect cannot be avoided. 
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